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CHAPTER I 
INTRODUCTION 

Spallation is a failure condition caused by stress 
wave interaction from blast or impact loading of an object. 
The most common example is the conical hole produced in 
glass by the impact of a pellet. In more ductile materials, 
such as aluminum, a spall may be as subtle as the formation 
of an internal void or as dramatic -as in glass, with the 
casting off of a chunk from the side of the plate opposite 
to the side impacted. 

Previous work on stress waves and spallation has 
been done analytically, experimentally, and numerically. 
There is a great deal of work on infinite half-spaces, but 
only a few studies on plates, which are applicable and 
necessary in order to understand spallation of plates. 

Davids [ 1 ] 1 has solved the problem of waves generated by a 
point load applied to a plate both before and after reflect- 
ing from the back face. A recent paper by Viswanathan and 
Biswas [2] studies waves in a plate produced by distributed 
loads. Their work is also analytical. Dally and Riley [ 3 ] 
have used dynamic photoelasticity to study stress waves in 

•'•Numbers in brackets, [], refer to the Bibliography. 
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a plate, but their work is not general. Ang and Newmark [4] 
studied stress waves by a finite difference technique but 
did not discuss reflections or spallation. Costantino [ 5 ] 
outlined the finite element method for studying stress 
waves but did not investigate the waves in detail or inves- 
tigate reflections and spallation. 

Figure 1 shows an experimentally produced spall in 
an aluminum plate. Obviously, large plastic and hydrodynam- 
ic forces were present in the total response of the plate to 
the impact. While a numerical scheme could be developed to 
approximate the total phenomenon, it is not necessary. The 
material must pass through the elastic range on the way to 
a spall condition. Investigation of the elastic waves 
alone will show which ones interact and where they interact 
to develop high stresses. These interactions and stresses 
may be interpreted for their contribution to spallation. 

The research presented here uses dynamic finite 
elements to investigate elastic stress waves in a plate. 

The investigation will discuss all waves produced by point 
and distributed loads applied to a plate. Conclusions will 
be drawn as to their importance in spallation. 
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c) Edge view 



d) Cross se.ction 

Figure 1. Experimentally Produced Spall in an 

Aluminum Plate 



CHAPTER II 


FINITE ELEMENT METHOD 

A short description of the finite element method in ' 
general is given first to outline the salient points of the 
method. Then a more detailed description of the specific 
method for axisymmetric problems is given. The discussion 
of the axisymmetric dynamic finite element method also out- 
lines the program used for this research. Detailed deriva- 
tions of the items mentioned are in Appendix A. The program 
is listed in Appendix C. 

General Method 

The general finite element approach to the analysis 
of any continuous body may be listed in seven steps. These 
are as follows: 

1. Divide the body into suitable pieces or 
elements ; 

2. Assume a displacement function for displacements 
in an element. Compatability is insured because 
the functions include parameters defined in 
terms of nodal displacements; 
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\ 

3. Nodes are the intersections of the imaginary 

\ 

lines and planes used to divide the continuum 
into elements. The geometric locations of the 
nodes must be recorded and also which nodes 
correspond to each element; 

4. Using force equilibrium at the nodes and virtual 
work, the elements are combined to a system of 
equations approximating the body. The equations 
involve stiffness, displacement, and applied 
forces; 

5 . If the problem, is dynamic, the inertia must be 
accounted for in the equations developed in the 
previous step. For transient wave problems, 
the mass is lumped at the nodes. Standing wave 
problems should be solved with the "consistent" 
mass matrix which is derived by virtual work 
with the inertia forces treated as distributed 
body forces. The mass matrix times accelera- 
tions of the nodes is added to the static 
equations of step’ (4); 

6 . The system of equations is solved for nodal dis- 
placements; 

7. Once the nodal displacements are known, they are 
used to find strains element by element which are 
in turn used to find stresses throughout the body 

In a static problem, the equations formulated in 
step (4) consist of stiffness terms times nodal displacements 
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set equal to forces externally applied to the nodes.. The 
nodal displacements are determined by an elimination tech- 
nique or by inversion of the stiffness matrix. In a dynamic 
problem, the mass times acceleration terms create second 
order, linear differential equations. These may be solved 
by integration, which may be prohibitively difficult, or by 
rewriting the acceleration in finite difference notation. 

The resulting algebraic equations are solved stepwise in 
time. The displacements, strains and stresses change with 
each time step. 

The lumped mass approach is based on a paper by 
Costantino [5]. The following axisymmetric application is 
based mainly on a discussion by Zienkiewicz [6]. 

Dynamic Finite Element Method for Transient 
Waves in an Axisymmetric Plate 

Element Formulation 

A circular, symmetric, distributed load applied to 
-an infinite plate produces an axisymmetric response in the 
plate. The infinite plate is approximated by an axisym- . 
metric disk of finite radius where the radius is great 
enough so that reflections from the radial boundary do not 
interfere with reflections from the back face of the plate 
until after the time of interest. This axisymmetric model 
is easily divided into axisymmetric finite elements as shown 
in Figure 2. The elements are triangular in the R-Z plane 
for simplicity of formulation and are the same as plane 
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strain elements with the exception that thickness of an 
element is replaced by centroidal circumference. Also, 
they must be capable of circumferential strain and stress. 
Nodal points of plane theory become nodal rings in 
axisymmetric application. Circumferential effects are a 
function of radial effects only, so that only radial and 
axial displacements are required to fully describe the prob- 
lem. Suitable displacement functions involve one constant 
for each degree of freedom of the element. This is a 
Ray liegh-Ritz approximation. 

UR = ot^ + + a^Z 

( 1 ) 

UZ = + a^R + a^Z 

where a.^ through are the unknown constants. 

R and Z are the coordinate position of the point 
whose displacement is desired. 

UR = radial displacement of the point. 

UZ = axial displacement of the point. 

The constants, o^, (i =1,2, . . . , 6), may be 
defined in terms of nodal displacements as is shown in 
Appendix A. The displacement of any point in an element 
is then defined by the nodal displacements and the geome- 
tric position of the point. 
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{e> e = [B] { 6 > e (3) 

where [B] Is a matrix of coefficients based solely on 

geometry of the element. The matrix {6} is the nodal dis- 

e 

i ■ 

placements of the element. 
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Hooke’s law for axisymmetric problems reduces to 
four stresses. 



Stiffness Matrix 


Element stiffness may be calculated by the princi- 
ple of virtual work. The matrix {S> e is the nodal forces 
resisting deformation of an element. Then from the princi- 
ple of virtual work: 

(S, e {a) I ’ fv 0EdV ' 

Substituting equations ( 3 ) and ( 4 ) and integrating, 

{S} o {6}^ = 2 it • RB AR •A , {6}^[B] T [D][B]{6} 

c c G C 

{S} = 2ir-RBAR-A-[B] T [D3[B]{6} . 

6 c 

RBAR is the radial distance to the centroid of the 
element. Therefore the stiffness matrix is 

[ K 3 e = 2Tt-RBAR-A-[B-] T [D][B]. , (5) 
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The stiffness for the system of elements represent- 
ing the continuous body may be obtained by adding the ele- 
ment stiffnesses in the proper way. This may be done in a 
tedious fashion by writing force equilibrium at each node 
of the entire system and also element by element. The 
terms of the element equilibrium equations are added into 
the appropriate positions in the set of equations for the 
entire system. This may be simplified by "globally" num- 
bering the nodes in any orderly sequence and then listing 
their displacements in a {6} array for the entire system 
sequentially from smallest node number to largest. Each 
element's nodes are called i, j, and m, numbering counter- 
clockwise (See Figure 3). Nodes i, j, and m have values 
prescribed by the global numbering. Then any number, k 12 , 
for example, in the element stiffness matrix adds into a 
position in the assembled stiffness matrix calculable from 
the node numbers i, j, or m which correspond to k^* 

The element stiffness matrix is 6 x 6 because the 
element has six degrees of freedom. The order of the 
assembled stiffness matrix is 2n x 2n where n is the num- 
ber of nodes and 2 is the number of degrees of freedom per 
node. However, the assembled stiffness matrix is banded 
and sparsely populated for the element system of this 
research. In fact, the maximum number of non-zero terms 
on any row is fourteen. This is because the greatest num- 
ber of nodes adjacent to any one node is six; so one node 
plus six adjacent times two degrees, of freedom is fourteen. 
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Great savings in core storage are achieved by making use of 
this fact. An additional matrix called NADJ, composed of 
the node numbers of each node and the nodes adjacent to it 
is used to condense the assembled stiffness matrix, KASY, 
from 2n by 2n to 2n by 14. A further explanation is given 
in Appendix A. 

Mass Matrix 

The lumped mass matrix is used in this research 
because of the assumption made by Costantino [5] that the 
consistent mass matrix would not allow wave propagation 
and in fact give rise to an infinite wave speed. That is, 
all points in the body would feel the disturbance at the 
same time. Appendix B shows that Costantino' s assumption 
is valid and compares the consistent mass matrix and lumped 
mass matrix in a simple closed form problem. The consistent 
mass matrix is derived in Appendix B. 

The lumped mass matrix has proven to give appro- 
priate results. The wave speeds are correct and stress 
distributions follow the correct shape and approximate the 
correct magnitude as is discussed in Chapter III. 

The lumped mass matrix does not require any deriva- 
tion. The mass of an element is equally divided among its 
three nodes. A node shared by six elements then has the 
mass of two elements concentrated at it. A node shared by 
three elements has the mass of one element and so forth. 

The lumped mass matrix is diagonal; Formulation of the 
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mass matrix for this axisymmetric problem is shown In 
Appendix A. 

As an example of the lumped mass matrix consider 
the mass times acceleration for the following one dimen- 
sional bar example. 



Equations of Motion 

Consider a simple spring and mass system. 



x 
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Statically; IF = 0 

-k6 + W = 0 
k 6 = W 


Dynamically ; 


IF = ma 


-kx + W = mx 
mx + kx = W 


This approach holds for a complicated system of 
masses and springs which is exactly the approximation of 
dynamic finite elements with mass lumped at the nodes. 

In matrix form, 

CM] {6 } + [K] { 6 } = {R} (6) 

where 



n is the number of nodes, and the matrix {6} is a columnar 
array of nodal acceleration arranged in the same order as 
{&}. 
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To solve this set of equations, the accelerations 
are rewritten in finite difference notation. Thus, 


U 


t >0 


U (t + At) ~ 2 U (t) + U (t - At) 

.2 


(7) 


( At )' 


Substituting the above relation into the equations 

of motion (6),- each equation may be solved algebraically 

for U/ , \ • The equations are elastically coupled, 

(t + At; 

but this does not interfere with the solution which is done 

node by node stepwise in time by increments of At . The 

time step At can be no larger than the time required for 

the dilatation wave to cross one element. The elements are 

DE 

right isosceles triangles so At ma x = *707 q^* > where DE 

is the length of one leg of the triangle and C-^ is the 

DE 

dilatation wave speed. For convenience At = .5 was 

used. This time step produces results indistinguishable 

from At max * 

The initial time step is a special case in which 
the condition 


U 1 


u 


t + At 


- u 


t - At 


t=0 


2At 


= 0 


( 8 ) 


must be used. 



CHAPTER III 
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ACCURACY OP PROGRAM 

Ang and Newmark [4] solved a problem involving the 
blast loading of earth with a finite difference technique. 
Costantino [ 5 ] used the same problem to check his finite 
element program which used rectangular elements. The pro- 
gram used in this research (Appendix C) used triangular 
elements (Figure 3) and produced similar results to the 
previous two solutions. Figure 4 shows the load distribu- 
tion on the surface. The similarity of the curves in 
Figure 5 shows each method is satisfactory. 

Results from the program in Appendix C were uti- 
lized to make Figure 6. The solid lines are theoretical 
wave positions calculated from the well-known formulas for 
wave speeds in a solid (see Mason [8]). 



Dilatation wave speed 
Shear wave speed 


(9) 

( 10 ) 


where X and y are the Lamd constants. 

Although the velocities computed are delayed in 
time, they are correct. The dilatation wave speed was 
located by the peak compressive axial stress generated by 



20 1 


18 



Figure 4. Cross-Sectional View of Grid and Load Distribution for 

Earth Problem 
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a sinusoidal point load, and the shear wave position was 
located by the peak radial displacements generated by the 
same sinusoidal point load. The delay in time is caused by 
the inertia of the lumped masses resisting the motion 
imparted by the forcing function. A sample print-out is 
included in Appendix C. 

Figure- 7 shows the temporal history of stress at a 
certain depth on the centerline. The theoretical solution 
is due to Davids [1]. An exponentially decaying step point 
load (Figure 8) was applied on the centerline, normal to the 
surface of the plate. The severe oscillations are due to 
the element size and the discontinuous nature of the forcing 
function. The discontinuity was reduced by applying a ramp- 
step load shown in Figure 9* Davids stated that for a step 
load the stress at a point would approach a value greater 
than zero for large times as shown in Figure 10. Figure 10 
shows; that the longer the ramp time, the smaller the oscil- 
lations. In fact it can be seen from Figure 5, that the 
oscillations become unimportant when the discontinuities in 
the forcing function are small. When a full sine forcing 
function is applied, the oscillations are unnoticeable as 
can be seen in Figure 11. The tensile peak at 1 micro- 
second is due to the shear wave. 

A more precise understanding of the nature of the 
oscillations may be achieved by writing the finite element 
equations for a one dimensional bar and solving them by 
Laplace transforms. This procedure is carried out in detail 



Axial Stress in KSI 
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Figure 7. Stress Response at a Point on the Centerline 
due to an Exponentially Decaying Point Load 



Figure 8. Exponentially Decaying Step Load 








Axial Stress in KSI 
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a. Full Sine Pulse 



Figure 11. Full Sine Pulse and Response 
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in Appendix B. The result is that at time zero for a bar 
approximated by two elements, using the lumped mass 
approach the reaction at one end is equal to the applied 
load at the other end. The same is true when the consistent 
mass matrix is used. This is because two elements are just 
one finite difference space for acceleration. For four 
lumped mass elements, the reaction is one-fifth the applied 
load. Using four consistent mass elements the reaction 
equals the applied load at time zero. The exact solution 
is that the reaction at time zero is zero. The more 
lumped mass elements that are used, the smaller the 
reaction. The consistent mass elements always give a 
reaction equal to the load, no matter how many elements 
are used, which indicates that the disturbance is propa- 
gated at an infinite speed. 

The reason for the non-zero (albeit small) reaction 
in the finite element approximation is that the solution is 
a Fourier series with a finite number of terms. A plot of 
such a series is a curve which oscillates about the exact 
answer. With enough terms in the series, the oscillations 
become indistinguishable, except at discontinuities. At a 
discontinuity, a jump known as Gibb’s phenomenon always 
occurs and is proportionate to the size of the discontinu- 
ity. The number of finite elements along the length is 
analagous to the number of terms in a Fourier series for 
approximation of waves which propagate axially. Therefore, 
the smaller the element size, the more terms in the series 
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and the better the approximation. Further, the first term 
in the series is the mean value of the series, hence the 
size of the reaction at time zero is just a measure of 
accuracy. 

As a final check of the program, a step load was 
applied to a long cylindrical bar. Figure 12 shows the 
results for two element sizes at about the same instant in 
^ time. The smaller elements give a more accurate answer in 
terms of percent error and rapidity of convergence behind 
the wave front. The small difference in the two results 
shows the larger element size is actually a good approxima- 
tion, and the results are influenced strongly by Gibb's 
phenomenon. The number of nodes in the two problems is 
equal. The only difference in the two bars is the element 
size. The grids are identical, each having five nodes 
radially and 168 nodes axially. The different element size 
changes the dimensions of the bar, but the number of terms 
in the approximation is the same. Therefore, the number of 
terms in the Fourier series for each element size is the 
same and the accuracy of each should be similar. The series 
for the smaller elements is somewhat more accurate because 
the coefficients of the terms in the series are better 
proportioned. 

It is clear from the previous assessment of accuracy 
that the lumped mass matrix does give satisfactory results 
with sufficiently small elements and continuous forcing 
functions. The accuracy improves with finer grids (more 
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computer space). While the accuracy of the stresses was 
limited by the capacity of the computer used, they do com- 
pare favorably with other numerical and exact solutions, 
and the necessary waves are generated at the proper speeds 
and locations to allow study of spallation. Limitations 
do exist, however, on the allowable sharpness of a dis- 
continuous input load. 
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CHAPTER IV 


RESULTS 


Although all four stresses (radial, circumferential, 
axial, and shear) are present, examination of the stress 
distribution has shown that the axial stress is dominant, 
especially in the area where a fracture may develop. 
Experimentally produced spalls are usually tensile failures. 
Since the axial stress is the most prominent, and undergoes 
the most dramatic fluctuations, and the applied stress is 
in the axial direction, the results shown are for axial 
stress. This will serve to illustrate the major response 
of the plate. It would be futile to attempt to convey all 
or most of the information generated by the program in the 
form of graphs. Each problem discussed produced approxi- 
mately forty stress distributions and nodal displacement 
displays which were each for separate instant in the 
response history of the plate. Condensation of these 
results to comments and a few graphs became essential. 
Direction of load application was always normal to the 
front face of the plate.' 
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Point Load 

The results of the point load application compare 
favorably with Davids [1] predictions. Although the step 
load response oscillates severely for reasons explained in 
Chapter III and Appendix B, Figure 6 shows the response 
does follow the proper values in time and space. The 
response is smooth when a full sine wave is applied as 
the forcing function (Figure 11). 

The full sine forcing function has a shape which 
produces a clear peak in the response which can be easily 
followed as the waves propagate. The lack of oscillation 
and the obvious peak made the full sine forcing function 
the -convenient choice for investigation of the stress 
waves. For the point load problem, a pulse of 100 pounds 
maximum amplitude was applied to the centerline in the 
axial direction. Figure 13 shows the axial stress dis- 
tribution for several times, and the significant features 
are labeled. 

The stress. on the centerline is the only one shown 
because the magnitude is highest on the centerline. 
Examination of the semi-graphic printout (see sample print- 
out in Appendix C) showed the stress waves moving out 
symmetrically from the point of load application in a pat- 
tern similar to the ripples produced by a stone dropped in 
a puddle. -Figure 14 shows this pattern in idealized form 
before and after reflection of the dilatation wave from the 
back face. Figure 13 shows the response of the plate 



Axial Stress in KSI 
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subjected to a full sine pulse whose period is about the 
same as the time required for the dilatation wave to cross 
the thickness of the plate. Figure 15 shows what happens 
when the period of the pulse is equal to the time required 
for the dilatation wave to cross five-eighths of the thick- 
ness of the plate. . 

Distributed Load Over a Small Radius 

A small radius is one which is equal to or less than 
about one-third of the plate thickness. The wave initially 
has a plane front, but seems to round off to a shape like 
that of the point load and produces a similar response. 

The full sine pulse was used for distributed loads, 
also. A maximum amplitude of 1,000 psi was chosen for 
convenience. 

The wave generated from the front face propagates in 
the axial direction for radii less than the radius of load- 
ing. From the point at the outer edge of the distributed 
load, the wave moves out like the waves produced by a point 
load. Since the dilatation wave moves away perpendicularly 
from the front face, it should not be followed by a shear 
wave. On the centerline it is not immediately followed by 
a shear wave. 

But at the extreme radius of loading the dilatation 
wave is followed by a shear wave immediately. ' It is produced 
because the dilatation wave has a radial component at that 
point. As time goes on, the shear wave spreads from the 
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Figure 14. Idealized Wave Patterns Produced by Point Load 


TTF = 2. 426 ysec 
E = 10.7 x 10 6 
v = .3125 
Y = .101 pci 
DE = 0.04 in. 



Figure 15. Axial Stress on Centerline Caused by Point Load 

TTF =2.426 

"A 
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radius of loading radially and axially. The radial propaga- 
tion is actually a radial dilatation wave caused by the 
Poisson effect. The reflection of this wave from the front 
face generates the shear wave. With the small radius of 
loading, the shear wave spreads to the centerline and fol- 
lows the axial dilatation wave. The smaller the radius of 
loading, the sooner the shear wave reaches the centerline. 

In addition, the part of the plane wave at the radius of 
loading attenuates causing the plane wave to apparently 
"round-off" to a shape similar to the shape caused by a 
point load. In comparison to a point load, the shear wave 
initiation on the centerline is delayed when the load is 
distributed. Figures 16, 17 , 18 and 19 show axial stresses 
on the centerline for four different radii of loading at 
different times. Again the centerline is the location of 
the maximum stress now because of the rounding off of the 
dilatation wave. 

Figure 20 compares the peak axial stress on the 
centerline as a function of depth for four radii of loading 
and a point load. The point load decays by the inverse of 
depth squared. The distributed loads approach a non- 
decaying plane wave as the radius of loading increases. 
Figure 21 shows the idealized wave patterns at three times. 

Distributed Load Over an Intermediate Radius 

An intermediate radius is, less than or equal to the 
plate thickness and greater than one-third of the plate 
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Figure 16. Axial Stress on Centerline Caused by Load 
over Small Radius, a = 0.08 in. 
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Figure 17. Axial Stress on Centerline Caused by Load 
Over Small Radius, a = .16 in. 
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Figure 19. Axial. Stress, on Centerline Caused by Load 
Over Small Radius, a = .32 in. 







thickness. The maximum axial stress develops on the cen 
terline, but high axial stresses also develop off the 
centerline. 

The high axial stress off the centerline becomes 
more pronounced as the radius of loading increases. If the 
axial stress distribution is thought of as a contoured sur- 
face, a saddle of lower axial stress exists between the on 
and off centerline regions of high axial stress. 

The on-centerline axial tensile stress develops 
first. The larger the radius of loading, the closer its 
location of development approaches a depth of about 70 per 
cent of the plate thickness. The" larger the radius of 
loading, the deeper the wave on the centerline travels at 
its applied stress value, resulting in less attenuation and 
higher peak tensile axial stresses for larger loading 
radii. The location of this stress buildup corresponds well 
with the location of th~ superposition of the edge and 
dilatation waves. This location is dependent on the wave . 
speeds which are functions of the material properties. 

The off-centerline tensile stress increases to a 
maximum located at a radius less than the loading radius 
but at a time well after the time the high tensile stress 
develops on the centerline. At the time the off centerline 
stress is maximum, it is the greatest tensile stress in the 
plate, but lower than the previously developed highest 
stress on the centerline. The location of the off- 
centerline maximum was observed to occur at half the 


ill 


loading radius for a loading radius equal to the plate 
thickness . 

Figure 22 shows idealized wave patterns for this 
case. Figures 23 and 2^ show axial stress distribution. 
Notice that the curves are for radii off the centerline as 
well as on it. The applied loads were uniform distributed 
loads. 

Distributed Load over a Large Radius 

When the loading radius is greater than the plate 
thickness, the off-centerline tensile axial stress develops 
first, propagates toward the centerline and exceeds the 
stress on the centerline. Once it reaches the centerline 
it continues to increase briefly, then attenuates. It 
stays on the centerline and propagates toward the front 
face. While developing and moving toward the centerline, 
the stress also moves somewhat toward the front face. 

After it reaches the centerline another "hot-spot" develops 
off the centerline with lower stresses between it and the 
centerline. The new hot-spot develops at the same depth as 
the current maximum centerline stress, moves with it toward 
the front face and also toward the centerline. This new, 
late developing hot-spot is smaller than the maximum center 
line stress. It is the same as the off centerline hot-spot 
developed by a load applied over an intermediate radius. 
Figure 25 shows idealized wave patterns for this case. 
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Figures 2 6 and 27 show axial stress at different radii and 
times. 

Table 1 shows important stress values, times and 
locations for each loading condition. Comparison of the 
times in the last two columns shows whether the maximum 
stress is due to superposition of the reflected dilatation 
(S-p with the incident shear (S^) or with the edge wave (S 2 ). 
The table summarizes several figures in this chapter. 
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l'M f :;uro 27. Axial Stress Due to Loading Over Large 

Radius, T = 5.27 psec 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 

Detailed, quantitative failure criterion are not 
discussed in this presentation, but a qualitative discussion 
is given on the nature of the response of an elastic plate 
to various impulsive loading conditions. 

Point Load 

The axial point load on the centerline sends a 
dilatation wave in all directions, but primarily (in terms 
of magnitude) in the axial direction. The Poisson effect 
causes the radial propagation which reflects from the front 
face as the wave expands. This reflection produces the 
shear wave, which follows the axial dilatation wave at the 
shear wave speed. The shear wave has particle motion trans- 
verse to the direction of propagation. The dilatation wave 
has particle motion in the same direction as propagation. 

The shear wave produces an axial stress opposite in 
sign from the axial stress caused by the dilatation wave. 
When the speedier dilatation wave reflects from the back 
face (a free surface), the reflection returns at the 
dilatation wave speed with axial stress of the opposite 
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sign. The reflection also produces a shear wave, but this 
shear wave occurs too late to influence spallation. The 
reflected dilatation wave and incident shear wave, there- 
fore, have axial stress of the same sign and superpose con- 
structively to cause a high tensile stress on the centerline 
This high stress can be greater than the tensile strength of 
the material and lead to a spall. Although the front face 
may be damaged by the impact or blast, this damage does not 
necessarily pierce the plate and, until the shear wave and 
reflected dilatation wave superpose, stresses are not great 
enough to fracture the material. For a point load, the 
crack initiates on the centerline and spreads radially. The 
location on the centerline may be found from graphs similar 
to Figure 6 for step loads and Dirac delta loading. How- 
ever, the sine pulse attenuates in a way such that the ten- 
sile stress develops at a shallower depth than predicted by 
Figure 6. 

Small, Uniform Loading Radius 
Initially, the stress waves differ from the point 
load, but after the shear wave develops on the centerline, 
the response is very similar to that of a point loaded 
plate. The dilatation wave has decayed away from the cen- 
terline so it seems to round off, but the front is a 
straight line radially. This attenuation also occurs on the 
centerline (after a certain depth) but to a lesser degree. 
The -stress will build up in the same way as for the point 
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load. The fracture condition develops on the centerline and 
spreads radially. Attenuation of the sine pulse results in 
the location of the tensile stress to be at a shallower 
depth than predicted by Figure 6. 

Intermediate Uniform Loading Radius 
The dilatation wave front decays inward from the 
radius of loading until it is "rounded" to the same condi- 
tion described for a small radius of loading. This occurs 
before the dilatation wave reflects from the back face and 
exhibits the same phenomenon which occurred in the waves 
generated by a uniform load applied over a small radius. 

Unlike the case of the small loading radius, a 
"hot-spot" of tensile stress develops off the centerline 
at a position and time indicating the cause is superposi- 
tion of incident shear and reflected dilatation waves. How- 
ever, this hot-spot develops late- and is smaller than the 
stress which develops on the centerline. 

The stress on the centerline results from the combi- 
nation of the reflected dilatation wave and the incident 
wave labeled S 2 in Figure 25 (see discussion for large load- 
ing radius). It develops on the centerline because the 
dilatation wave has attenuated away from the centerline. 
Therefore, the fracture develops on the centerline and prop- 
agates outward. A secondary, crack may develop off the cen- 
terline if the hot-spot has a high enough tensile value. 
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Large Uniform Loading Radius 
The numerical results from a large loading radius 
clearly support the mathematical analysis of Viswanathan 
and Biswas [2]. The highest tensile stress develops off 
the centerline before a high tensile stress develops on the 
centerline. Also, 'it develops too early to be caused by 
combination of the shear wave and reflected dilatation 
wave. It moves toward the centerline while moving toward 
the front face. It does not develop until after reflection 
of the dilatation wave. These conditions mean that • the 
reflected dilatation wave is combining with what Viswanathan 
and Biswas call "an edge wave . . . with a toroidal front, 
expanding with velocity c," the dilatation wave speed. 
Figure 25 shows this "edge wave" labeled "S 2 ." "Edge" 
refers to the outer edge of the loaded area on the front 
face, from which the edge wave emanates. Therefore the' 
crack initiates off the centerline and propagates toward it. 
The location- of the off-centerline tensile stress develop- 
ment is at a -radius less than the radius of loading due to 
/ 

the same attenuation "round-off" described earlier. 

The stress on the centerline may exceed the tensile 
strength of the material before the off-centerline crack 
gets there. A second crack would develop on the centerline 
and travel outward to meet the first crack in this case. 

Finally, it should be noted that attenuation of the 
wave resulted in the four cases described. If the wave 
front had maintained the input stress between the centerline 



and radius of loading, only the two cases of point load and 
distributed load need have been discussed. Without atten- 
uation even the smallest of loading radii would produce the 
highest stress off the centerline first due to superposi- 
tion of and S 2 - This would occur near to the back face 
and at the radius of loading. The shape of the waves does 
not attenuate or round-off. The magnitude does attenuate 
from the radius of loading (least attenuation on the center 
line), making the waves seem to round-off.; The reason for 
this attenuation may be thought of as a "diffusion" of 
stress from the highly stressed (saturated) region inside 
the radius of loading to the region outside the radius of 
loading where the plate has low stresses. This diffusive 
process is due to the Poisson effect which causes wave 
propagation in directions other than the direction of 
loading. 
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APPENDIX A 


DERIVATION OP FINITE ELEMENT MATRICES 


Displacement Functions 

The displacements within elements are given by: 


UR = ot-^ + oigR + a^Z 
UZ = ajj + a^R + agZ 


(A. 1) 


where R and Z are the coordinates of the point whose 
displacement is being calculated. UR and UZ are the 
radial and axial displacements, respectively. Since the 
displacement functions are valid for the nodes as well as 
inside the element, the unknown coefficients through 

ag may be determined by writing the displacement functions 
at each node. Let the nodes and elements of the system be 
numbered in the manner of Figure 3. Let the nodes of each 
element be i , j and m in a counter-clockwise manner and 
refer to the global numbers of Figure 3- 
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UR ± 

= 

a l 

+ 

°2 R i 

+ 

“3 Z i 

UR. 

1 

= 

a l 

+ 

“2 R j 

+ 

“3 Z J 

UR 

m 

= 

“1 

+ 

“2 R m 

+ 

“3 Z m 

•H 

!SJ 

= 

a i| 

+ 

“5 R i 

+ 

a 6 Z i 

UZj 

= 

a 4 

+ 

“5 R j 

+ 

6 3 

UZ m 

= 

a 4 

+ 

a c R 
5 m 

+ 

a/-Z„ 
b m 

The 

constants 

a l 

> ct 2 


(A. 2) 


the first three equations by Cramer's rule and , a,- , 

and ag are found in the same way from the last three 
equations. One example and other results follow. 


“i = 


ORi Z ± 

UR J R J Z J 


UR R Z 
m m m 


i Rl Zl 


1 R J Z 3 


1 R Z 
m m 


(A. 3) 


Solving A. 3 gives: 

(a. UR, + a. UR, + a. UR m ) 
i i .1 .1 mm 

a-, = rr ; — t , — r — s 


(a i + a j + V 


(A. 4) 


where (a. + a. + a )• = 2A = twice the area of the triangle 
i 0 m 

Similarly the other a's are defined in terms of nodal 
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TABLE 2 


COEFFICIENTS USED IN THE DEFINITION OF DISPLACEMENT 



The definition of strain in terms of displacement 
for axisymmetric problems is:' 



II 

K 

u 

8UR _ 

8R °2 

60 


UR “l , . C, 3 ZBAR 


£ 0 = 

R " RBAR a 2 RBAR 


e z = 

8UZ 

8Z " a 6 

(A. 6) 

y rz 

8UR . 3UZ 

3Z + 3R °3 + “5 



The definition of e Q includes the following approx 
imation since these are constant strain elements. 


RBAR = 


ZBAR = 


< R 1 * R ,j + R m> 


< Z i + Z 1 + V 


(A. 7) 


Define the element displacement matrix as: 

UR j 
UZ, 


{6} e H 


UR 


UZ. 


UR_ 


m 


UZ 


m 


(A. 8) 


Substitution of the definitions of the constants, 
a^, results in the strain displacement equations in matrix 
form. Equations A. 6 become: 
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f e 1 
e R 


b i 

0 

b J 

0 

b 

m 

0 

£ 0 

1 

d i 

0 


0 

d 

m 

0 

C Z 

2A 

0 

c i 

0 

c j 

0 

c m 

r 

5d 

1S1 


C i 

b i 


b J 

c 

m 

b m 


where : 

c^ZBAR 

d i = RBAR + b i + RBAR " 

a. c.ZBAR 

d j = RBAR + b j + RBAR 

a e ZBAR 

^ _ m j_ v, j_ _m 

d m " RBAR D m RBAR 


(A. 10) 


The Stiffness Matrix 

The element stiffness was derived on page 10 to be: 


[K] = 2 • it • RBAR* A[B] T [D] [B] (A. 11) 

e 

In order to develop a procedure for assembling the 
stiffness matrix for many elements, consider a two element 

example, shown in Figure 28. 

For the system, equilibrium of nodal forces gives: 
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S 1 (1R) = 0 

S 1 (2Z) = R 1Z 

S 1 (2R) + S ? (2R) = 0 

S 1 (2Z) _ S 2 (2Z) = R 2Z 
S 2 (3R) = 0 

S 2 ( 3 Z ) = ^32 

s 2 (4r) + s 2 (4r) = 0 

S 2 (4Z) + S 2 (4Z) = Q ijZ 


Notation 



(1R) 

*/ 

Direction 
Global node number 
Element number 


(A .12 ) 


The element stiffness equilibrium equations define 
the nodal forces (S) in terms of stiffness times displace- 
ment. Substitution of these terms in the above system 
equilibrium equations results in the assembled stiffness 
matrix for the two element"' system. In order to show that 
the position of the element stiffness term in the assembled 
stiffness matrix is given by node number, the element stiff- 
ness matrices for the two element example are given in the 
equations below in expanded form.. The program in Appendix B 
uses the same minimum space (KELM) for each element stiff- 
ness matrix and in the subroutine KFORM computes the posi- 
tions to add them to in the assembled stiffness matrix 
(KASY) by node number and a pointer matrix (NADJ ) . Maximum 
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storage compression is achieved by using NADJ to eliminate 
the zero terms in a row, resulting in KASY having only 
fourteen columns. 


[K] 1 (5} 


CK] 1 = 


[K] 2 = 


1 — 1 

1 1 

ru 

{<$} = {R} 



k (11) 

K n 

k (12) 

k 12 

0 

k (1 4) 

13 

k ( 2D 

21 

k ( 22) 

22 

0 

k (2 4) 

23 

0 

0 

0 

0 

k (4D 

31 

k (42) 

32 

0 

k (W 

k 33 

0 

0 

0 

0 

0 

k (22) 

11 

k (2 3) 

12 

k ( 24) 

k 13 

0 

k (32) 

K 21 

k (33) 

23 

k ( 34) 

k 24 

0 

k (42) 

K 31 

k (i!3) 

32 

k ( ^) 
33 J 


(A. 13) 


(A. 14) 


(A. 15) 



Each k in the above expanded element stiffness 
matrices (A. 14 and A. 15) is a 2 by 2 partition of the 
original 6 by 6 element stiffness matrix. The subscript 
refers to the row and column position of the partition in 
the 6 by 6 element stiffness matrix. The superscript is 
determined by the i, j, m node numbers of the element and 
is the position of the partition in the assembled stiffness 
matrix. This is further illustrated below. Consider any 
one element with some external loading (R). 



(A. 17) 


(A. 18) 
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where , Fh , etc., are 2 by 1 partitions such that: 

f UR i) 

6^ = 1 V and so forth. 

Kj 

The system of forcing KELM to conform to the i, j, 
m sequence of nodal displacements in the {6} array 
insures that the proper indexing of the partition will occur. 

The Lumped Mass Matrix 

The lumped mass matrix is computed by the subroutine 
MFORM. Since the grid generated by this program is regular 
and defined by the number of nodes axially and radially, 
the mass matrix is formed node by node. The nature of the 
grid (Figure 3) allows the matrix to be computed in stages, 
as indicated by the comment cards (see Appendix C). For 
example, all interior nodes are shared by six elements, so 
the factor CF is equal to 2. The average of the radii to 
the centroids (RBAR’s) of the six elements surrounding any 
interior node is always the same as the radius of the node, 
so DR = 0 . Similarly the mass lumped at each node is 
based on the number of elements which share the node and the 
average of the radii to the centroids of those elements. 

The lumped mass matrix for the two element example in 
Figure 28 is: 



where , 


{M] = 


• 333pAC 1 0 0 0 0 0 0 0 

0 .333pAc 1 0 0 0 0 0 0 

0 0 .667pAc 2 0 0 0 0 0 

0 0 0 .667pAc 2 0 0 0 0 

0 0 0 0 .333pAc 3 0 0 0 

0 0 0 00 .333pAc 3 0 0 

0 0 0 0 0 0 .667pAc 2 0 

0 0 0 0 0 0 0 .667pAc 2 


c 1 = RBAR for element 1 

c 2 = R 1Z “ R 2Z 
c ^ = RBAR for element 2. 
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(A. 12) 
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Pa 1/67$,, 


vr 


t* 
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APPENDIX B 

DERIVATION OP LUMPED AND CONSISTENT MASS 
MATRICES FOR ONE DIMENSIONAL PROBLEM 
WITH LAPLACE TRANSFORM SOLUTION 

Description of Problem 

In order to determine which mass matrix approach 
should be used for propagation problems, both will be used 
to generate equations approximating a bar subjected to a 
step tensile load at the free end and fixed at the other 
end. Laplace transformation will be the method of solution 
of the equations resulting in an expression for the fixed 
end reaction at t = 0(+) • One dimensional finite ele- 
ments are used to approximate the bar which is assumed to 
have a length many times greater than the diameter. 

Lumped Mass Matrix 
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The subscript refers to the node and the superscript 
refers to the element. Let an element be described in 
general as follows and choose a linear displacement 
function. 

A = area 

E = Young's modulus 
p = mass density 



u = + a^x 


_ u .i U i 

e x a 2 


[B] = 4 [-1 1] 


[D] = [E] 



= 



C K ] e = A &[Bj[D][B] 


CK] e = f 


1 -1 

-1 1 


{6 >e = [K] e {6} e 


(B.l) 

(B. 2 ) 

(B. 3 ) 
(B.H) 

(B. 5 ) 

(B. 6 ) 
(B.7) 


(B.8) 
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°1 

1 2 

Sg + 


= R ] 
= 0 

= Q 


(B.9) 


AE 

i 


" 1-1 cf 


u l 

. 

R 1 

-1 ■ 2 -1 

< 

U 2 

s= * 

0 

0 -1 1 


u 3 
l -i 


Q 

L J 


[M] = m. 


0 

0 


I 0 


1 

0 


m = pA£ = mass of 
0 one 

element 


(B.10) 


(B.ll) 


[M] { 6 } + [K] { 6 } = {R} 


(B . 12 ) 


Let Q ( t ) = QH(t) . 


By following the pattern shown for two elements, 


the equations for four elements may be written: 
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o 


1 

2 

0 

0 

0 

0 


“l 



1 

-1 

0 

0 

cf 


U 1 

0 

1 

0 

0 

0 


U 2 



-1 

2 

-1 

0 

0 


u 2 

0 

0 

1 

0 

0 

< 


>+ 

AE 

l 

0 

-1 

2 

-1 

0 

4 

U 3 

0 

0 

0 

1 

0 





0 

0 

-1 

2 

-1 


u 4 

0 

0 

0 

0 

1 

2 


*5_ 



0 

0 

0 

-1 

1 

J 


u 5 

J 


R 1 

0 

0 

O' 

Q(t ) 


ra 


(B. 13) 


Divide both sides by , recall = pAl , and 


longitudinal wave speed, c = /E/p . 


“i + ( 

!) 2(2u i - 

2u 2 ) = 

"“1 

m o 

2u 2 + 

(!) 2 < 

- 2u l 

+ 4u 2 

- 2u^ ) = 

2u 3 + 

r< 

-2u 2 

+ 

- 2u ft ) = 

2u 4 + 

r< 

-2u 3 

+ 4u 4 

- 2u 5 ) = 

V ( 

l) 2 <- 

2u 4 

+ 2u 5 ) 

2QH ( t ) 


(B.14) 


Initially, at t = 0 , the displacements and velo- 
cities of all five nodes are zero. The boundary condition 
exists such that the left end is fixed, or, the displacement, 
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velocity and acceleration of node one are zero for all 
time. This boundary condition means that the Laplace 

• 4 

transforms of and u^ , are also zero. The Laplace 
transforms of the equations with the initial and boundary 


conditions substituted are: 


-(f) 

2 v 

L[R X ] 

m o 

s 2 U 2 

+ ( 

l) 2(2D 2 - V 

s 2 U 3 

+ ( 

l) 2 (-U 2 + 2U„ 

s 2 U 4 

+ ( 

l) 2( -°3 + 2U H 

s 2 U, 

+ ( 

h + 2U 


= 0 

- U 4 ) = 0 (B.15) 

- u 5 ) = 0 



where the capital U indicates the Laplace transform of 
the displacement and s is the Laplace operator. 

It is clear from the first equation that U 0 can 

*! C. 

be expressed as a constant times the transform of . 

This is substituted for U 2 in the succeeding equations. 
Then, with that substitution in the second equation, 

is expressed in terms of the transform of . The 

substitutions continue and after substitution and rear- 
rangement the fourth equation expresses in terms of 

the transform of R^ . The fifth equation is used to find 
the transform of R^ in terms of the transform of the 
forcing function. The result is separated by partial frac- 
tions into recognizable transforms which are the transforms 
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of sine functions and the transform of the forcing function 


all multiplied by constants. The inverse of this Laplace 
expression is the reaction, R^ , as a function of time. 
At time zero, the four sine terms vanish, and 


t=+0 


Q 

5 * 


(B.16) 


Additional elements result in more terms in the 
series and at time zero approaches zero. 

Consistent Mass Matrix 

The mass of the element is considered to be dis- 
tributed. Its inertia is assumed to act as a body force 

causing nodal forces, S b , in combination with the nodal 
0 

forces, S , caused by the stiffness of the element. The 

virtual work expression must then include the contribution 

of the mass intertia. 

.* 

Let X = pu = body force. 


S'lu. + S^u. 
i i 0 J 



ceAdx + 



XuAdx 


(B. 17) 


The displacement function is now variable with position and 
time. 


u ( x , t ) = a 1 + a 2 x(t) (B.18) 



(B.19) 



u(x,t ) 




rewriting B.17 


(S b + S e )J u.^ + (S b + S e )j Uj 


rXj 


aeAdx + 


XuAdx 
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(B . 20 ) 


(B . 21 ) 


The contribution of the inertia forces may be formed 
by separating the body force terms. 


<S b )i u, + (S *) 1 Uj 


1 

& 


r x 


j XA[Uj (x - x.^) - u i (x - Xj)]dx 


(B . 22 ) 


Equating coefficient of u^ and u. ; 


( sb )i = 7 


< sb) 5 - t 


fX 


^ XA(Xj - x)dx 


x 


x . 
J 


XA(x - x^)dx 


(B. 23) 


Substitution of X = pu, 

x ± = (1 - 1)Z , 

Xj = it. 

rearranging, substituting Z = 
ponding changes in the limits of integration, the nodal 
forces due to inertia forces may be written; 




and with the corres- 
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(S b )j = -m o 


-1 


( ( u j - u i )Z" ; + Uj Z )dZ 


(B.24) 


(S 6 )] = -(S b )J - m D 


-1 


( (Uj - U 1 ^ Z + U j 


Evaluation of the integrals yields; 


<S b )j 


(S b f 

J 


= m 


= m 


u . u. 

^ + T 


O 


J. 

3 


u . 

^ + TT 


(B. 25 ) 


By writing equilibrium of forces at the nodes of 
the two element problem, remembering that the nodal forces 
consist of both inertia and elastic forces, the following 
set of equations is formulated. 



— - 


* 

•• 


1 1 n 

3 5 0 


u l 


12 1 


u Q 

m 

o 

F 3 F 


d 1 


nil 


U-J 


0 F 3 


3 


+ 


AE 

l 


" 1 -i o’ 


U 1 


r 

R 1 

-1 2 -1 

• 

u 2 

ii 

0 

0-11 


u 3 


Q _ 


(B. 26) 


where the first matrix is the so-called consistent mass 
matrix and the stiffness matrix is the same as before. 
M6 plus K<5 give the nodal forces. 
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After formulating the four element problem with 
the consistent mass matrix, the solution of the five result- 
ing equations by Laplace transform follows the same proce- 
dure as shown for the four element lumped mass problem. 

The same initial and boundary conditions are used and the 
reaction at node one is found by successive substitution. 
Again, the solution is a step plus sine functions all multi- 
plied by constants. At time zero, only the step acts on 
node one and the constant coefficient is unity. The 
reaction at t = 0(+) is: 

= -Q (B.27) 

t=+0 

Even with additional elements, R 1 = -Q at time 
zero. Therefore, the consistent mass matrix does not allow 
propagation of transient phenomenon. The reaction, R^ , 
is a Fourier sine series in which the first term is a con- 
stant since the forcing function is a unit step. It is 
well known that Fourier series solutions for a step, or any 
discontinuous function, have a jump at the discontinuity 
known as Gibbs phenomenon and proportional to the size of the 
discontinuity (See Biot and von Kdrm&n [9]). Further, the 
above analyses indicate that the dynamic finite elements 
solutions are equivalent to Fourier series solutions even 
though the process is numerical in the computer. There- 
fore, the oscillations in the finite element solution are 



>- 
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due to too few terms in the series (too few nodes) and ( 
Gibbs phenomenon. 

Finally, it is easy to see that the equations of 
motion developed by lumping the mass at the nodes are the 
same as the one dimensional wave equation. The comparison 
is obvious when the - accelerations in the finite element 
method are written in finite difference notation. Then 
the differential equation for one-dimensional waves is 
written in finite difference notation. The finite element 
stiffness times displacement terms are already equivalent 
to the difference notation. However, when the accelerations 
in the finite element equations for the consistent mass 
matrix are written in difference form, the result does not 
correspond to the difference form of the one dimensional 
wave equation. This qualitative look at the nature of the 
equations further supports the idea that the consistent 
mass matrix will not allow propagation. 



APPENDIX C 
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APPENDIX C 

COMPUTER PROGRAM JMMSPALL 

Following is a listing of the program used in this 
research. A sample printout is included, consisting of a 
typical printout for one time step. In use, the program 
solves a dynamic axisymmetric problem and prints out 
stresses and displacements at several times. The sample 
problem is a plate with a distributed load. There are 
twenty-four nodes axially and thirty-five nodes radially. 
Distributed loads are produced by resolving the desired 
stress to nodal forces. 
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0.C9 0.09 3.1? 0.1 1* 3,1* 2.16 0.19 3.23 0.28 1.3* 0,*2 9.31 C. 6 l t , 69 Cn7t 9.82 ft. ©5 3»8* 0.79 3,7', 7.39 2-*7 1* 

-0.01-0. 03— 0.C6— 9.0 7— 3, C 8 -3. J 7 -0. ?5— 9, 01 O.C* 0.19 0»2e P.-.5 0.63 “.61 *>.*7 *,12 l.i> 1.27 1.12 l .“8 C,»b '.5* 3 . 16 

0.10 0.10 n.ll C.12 3.1* P.17 0.19 9.23 0.28 3.3* ft.*? 3.92 0 , 6 l C,7C 0.76. P.32 P .86 9,65 0 . 6 J 0,71 0,59 0 *6 2*35 
-O.Cl“0.0*-0.06-0.9B-3. 08-3. 37-3. 0*->.PO 0.C6 C.l* 0.28 C .*6 U. 6 * 0.62 <',98 1.13 1.26 1.29 1.2* 1,11 ?» 6 B 3,9* 3,13 

ICfNT" 0.02 0.06 9.10 O.l* 0.16 3.22 0.26 0*30 3.3* 0.36 0.*2 0.*6 0.50 p. 9* 0.56 0.62 9**6 0.70 0.7* 0.7b 0.82 3.66 2,90 


jot 


EACH BLOCK is A STPCS5 I* KSI 
2 STRESS IS KSI 


N004L 015®lttLKFNTS 


t« o.*632?&e-os 



HC6HT 

1.36 

1*32 

1.26 

1.14 

i.2© 

1.16 

1.12 

1.08 

1.06 

1.00 

0.96 

0.92 

0.80 

0.64 

0. 00 

0.76 

0.72 

0.68 

0.64 

0.60 

0.56 

0.52 

0.46 

0.44 

0.40 
0.36 
0 .12 
0.28 
0.24 
0.20 
0.16 
0.12 
0.06 
0.04 
0.0 


0.0 —0*8 -0.3 

0.6 

1.8 

3.3 

5.1 

6.0 5.4 4.9 

4.5 

4.1 

3.7 

3.5 

-1.1 -l.P -0.5 

•0.5 

1.7 

3.) 

5,2 

5.9 5.7 6.6 

5.3 

3.4 

5.: 

5.1 

-1*7 -1.2 -0.5 

0.5 

1.8 

3,4 

5.3 

0.5 6.4 6.5 

6.5 

6.6 

6,6 

6.6 

-2.2 -1.3 -0.4 

0.7 

2.9 

3.3 

3.4 

7.5 7.4 7.5 

7.6 

7.8 

7.9 

e.* 

-2c5 -1«* -0.3 

0. B 

2.2 

3.0 

5.5 

6.6 6.7 8.8 

8. 9 

5.1 

9.2 

9.3 

-2.6 -1.4 -0.2 

l.C 

2.4 

3,9 

5,5 

10. 0 10.1 10.2 

10.3 

13.5 

10.5 

1 J.5 

-2.5 -1.3 -C.l 

1.1 

2.4 

3.9 

5.* 

11.4 11*5 11.6 

11.7 

11.7 

11.7 

11.5 

-2.3 -1.2 -0.1 

1.1 

2.4 

3,7 

5,1 

12.7 12.8 12.9 

12.9 

12.8 

12.5 

12.4 

-2.C -1.1 -0.1 

1.0 

2.2 

3.5 

4,7 

13.7 13.0 13.8 

13.7 

13.6 

.*.3.3 

13.0 

—1.8 -1.0 -0.1 

0.9 

2.0 

3.1 

4,3 

14.4 14.5 14.4 

14.3 

14.1 

13.9 

13.5 

-1.7 -1.1 -0.3 

0.7 

1.7 

2.7 

3*8 

15#C 14.9 14,9 

14.7 

*.4.5 

14.3 

13.9 

-1.8 -1.2 *0.5 

0.4 

1.4 

2.3 

3s 3 

13.4 15.3 15.2 

15,1 

15.0 

14.7 

14.5 

-2.1 - 1.4 -:,7 

0.2 

1.1 

2.0 

2*9 

i:.9 15.8 15.7 

15.6 

15.5 

15.3 

15*1 

-2.4 -1.7 -1.0 

-J.l 

0.8 

1.7 

2,6 

16.4 16.3 16.2 

16.2 

16.1 

16.3 

15*9 


\ **.v *-*« ' - 

5.1 5.1 5.1 5 ,2 5v5 5.8 7 

7.2 9.? fl..O >2.* '3 9 14.7 .4*9 

6.6 6.6 6-3 6.4 6 5 6,. 3 6..* «•» 6 

7.3 9.1 »P-9 12.4 13 6 14.4 ?.6«- 6 14,6 !J.o . C,* 

P.U 7.9 7-7 2.4 7. I 6.9 6.7 6.e 6*9 7,2 


7.2 6.9 li.9 13.1 13.9 '.4.2 14.: 13.4 U,3 

4. 2 9.0 J* 7 8,3 0 7.4 7.i 7.4 T - ’ ** 


11.2 1C.B 10.3 9,7 9 1 8.6 8<*» h » * 

6.4 7.6 6.8 9* 9 19- B 11.7 »2*3 52.5 12.3 

12. J 11.5 10. 9 l ie 3 9.7 9.3 9. i 9., 9,6 9,9 

5o 9 7.0 8,1 9.1 1C. 1 11.0 U.7 12.2 12,7 ll»5 

12.6 12.0 11.5 tv. 9 Itt.S 10.2 iC.,1 l'«3 10,7 .4,3 

5.4 6.4 7*4 8« 3 9*3 10.3 11.1 11*6 11,7 31,4 

13.1 12.6 12-1 11.6 11.3 it*2 11*8 11*6 12,2 -3, u 

* o «-• K. 7 7,7 6.7 9.6 19.7 Mtl li>» '1,2 


13,2 

8,0 

*;.o 

0,7 

10,6 

9,1 

9« i 
9.2 

17,6 

7,3 

12.3 

7,8 

10, 8 

0, 1 

9 3 
8,3 

t 3.7 
7. ** 

1?,3 

7,3 

19.9 

7.6 

9.3 

7,b 

! 3.6 
6*9 

12,4 

7,2 

n,o 

7.5 

9,4 

7,9 

J 3,4 
7.2 

W,3 

7,4 

11.0 

7.6 

9,5 
"rJ ' 

» x , * 
7,7 

12* • 
0, - 

12,9 

9,5 

a <f 

12,7 

8,4 

11,9 

3,8 

5 

9,4 

9,7 

1 

12,3 

’ 1,7 

10,8 

9w 7 


5 0 6-9 

6,9 8,1 


6 7 3.-0 

4.6 C.9 


g.J 7.. 7 6 8 i. I 6, 3 

7.9 6,. 7 56 6.4 6.1 
7, 9 7 S 7. t 7 i 4 

7..:> 45 * 7 7,4 

6.1 » v 4 B 7 9. 2 9. 7 

6-3 7.5 7 j 7,7 6.6 

6 » 7 9,2 9 7 - ■*. ; 1 I C- • 6 


j * 7 9,8 

j* :U? lli 7 


11,1 

15,7 


10,7 9,9 

>2,2 &3.i 

1 0, 6 9,9 

13,9 15c) 

1 C a 6 9.9 

15.6 17,0 

10.5 9*9 

17.6 19*0 

10,5 9„ B 

i 3 9,7 2' -9 

i 13,3 9,i 

i 21.- 7 22-9 

10,2 9*5 

2 3,5 24.9 


H .6 

9 .6 

?, <j o- 4 0,7 9,5 11.0 

,0.8 1 .« 6 ’2 iZ<l 13.-0 

6- j $,8 9 2 iJ- 2 11 .7 

.2,3 13.1 *3:3 -4.3 »4; 6 

9.2 9.1 9 7 12,8 12,4 

j i. • 1?*1 15.7 *6,2 16 >4 

9-4 9* .2 : llifc 12 o 8 

it:. - Uc* l T 7 IS.* 18*3 

9 4 9,5 *' i !. »4 13,*. 

;p ; . :9,3 *9, 7 ::-i 2:*2 

9-» VS 3 lie 6 l!:3 
?.c • -8 27*1 

9.4 9.5 *•? 3 U.6 1? » 3 

22,1 2 jo 1 23 9 24.1 2‘.1 

9.2 9,4 J ’ 2 11.5 13,2 

, : 4 . «■>, i 2* 8 *6 * 26* 0 


... -0,* ;,J i,» 2. * 3.3 *.2 5.1 6.* 7.7 ».t ’-.l i|>i '»•’ 

17H ItIo 16*9 16*9 1 6o 8 16.0 16. 3 16.7 l6„7 16. « 17.1 .7.6 16.3 9*1 '* - * *' 

. , , . . . . . a . • , 1.2 3.1 4,?. 5.2' 6,5 7.- 8 9.1 IC.k li’t b 11* ) '“SJ 

In» wle ItIo 1 7. T 17.7 17. t 17-. e 17.8 17.7 19.1 19.5 17 1-17.7 r-.f i: .1 22.7 .7. 

in*2 iP.8 10*9 ?->»5 9,7 9, 

22.5 I3o5 24,6 25*7 27,0 2o, 


9 C 3 

17. : t . 


:ici 12*0 

18,1 27,9 


3 V > 0 9*3 

?5: 3 2eoi 


-4.1 -3.1 *2.'» -UP 0.1 i»l 
19.6 i*#.b 19. e 19,6 19.9 20#1 


2,1 

3.1 

4.1 

5.3 

6.6 

8.0 

15-9 

19.0 

19.2 

19.3 

20. 0 

20 >7 

2,1 

3a 1 

4.3 

5,5 

6. 0 

8. 2 

20.2 

20.4 

Z< .7 

?!.l 

?’c5 

?:*3 

2.2 

3.3 

4.5 

5,8 

7.2 

8*. 5 

21®7 

22.0 

22.3 

77, R 

23.3 

24,. C 

2.5 

3.7 

4.9 

6,2 

7.5 

e. 7 

23.4 

23.8 

24.1 

24,6 

25.0 

25 6 

2,8 

4.0 

5.3 

6* 6 

7,8 

8.8 

25.4 

26.3 

26.1 

26.4 

26.7 

27a 2 

3.2 

4.4 

5,6 

6.8 

7.9 

8. 8 

27,6 

27.8 

26.9 

20.2 

28c 3 

28.6 

3c 6 

4.T 

5.9 

6-9 

7,9 

8 6 

29,8 

29.0 

29 v 8 

29.7 

29.7 

2°c 7 

3.8 

- 4.9 

5.9 

6.9 

1 « • 

6. 3 

31,7 

31.5 

31-3 

3C® 9 

30.6 

3C 4 

3.9 

4.9 

5.8 

6. 6 

7.2 

7.7 

33.2 

32.8 

32.3 

31.7 

31.2 

30o 9 

3.8 

4.7 

5.5 

6.1 

6.6 

7.:0 

34.0 

33*4 

32.8 

32.1 

31.4 

31 *0 


9.4 ’0.3 i-%e i:.3 : >.3 


9,4 3*7 

26,6 3 . 


10.3 if- 7 

10,5 

9,9 

9,1 

8,3 

:5c 7 26,o 

27,6 

20,8 

30,1 

31.5 

*.Oa 3 13c 5 

;:,2 

9,6 

! R,7 

7 4 

27.. 1 28* j 

25,9 

?:,i 

31,4 

32,9 

1*a1 1 ■**> ? 

9,8 

9,1 

8,3 

7, 5 

2b»4 29*2 

30,1 

31,2 

32,6 

34 0 2 

9 V F 9.7 

9,2 

3,6 

7*8 

7.1 

29*4 3*«1 

31,2 

32.2 

33.7 

35,4 

9,3 9c 2 

8,7 

8,1 

* 7,3 

6,6 

3C.2 30.9 

34,8 

33. 3 

24,6 

36.4 

8a 7 2.5 

8.1 

7,5 

6,6 

6s> 

*>.0 31*. 4 

32,4 

33. T 

35s 4 

37 c 4 

8cf. T c 0 

7,5 

7».< 

6,3 

5v7 

31a* 31.6 

32,8 

24,3 

36.1 

38.2 

7-2 7.1 

6,8 

6.4 

5.8 

5c 2 

31.2 32.0 

33,2 

54,8 

26,6 

38.9 

6,4 6.3 

6,1 

5,6 

f , 3 

4.8 


3:2 ,K ,E ,K £2 £? 3?:! ,1:1 ,?;l a: 3!:*, 3*3:1 ,2:2 

3:2 &5 ,2:2 3 l 5 :i ,u: ,2:2 >2:5 ,?:! ,t :l «:i 30:: 3?:; 3!:$ 3^ *i:I 

,S:i ,?: 1 ,1:’, ,i:: ,i:i .n: J:: 3 ^.i ,5:: »:! 

0.6 0.6 0.8 1.3 1.7 -2.1 2.5 


6.9 9-, : •? p 12.-8 

27.3 .*•‘,9 *.9 a 50 . v 29-. 6 

6,6 3*7 9 5 3*8 12.4 

29:6 3 i. 7 J 1 . 5 ?» -8 31,6 

5 2 9 . 13.4 12- > 

Jl,,i J2..4 SV 3 !-3 6 25- 3 

7.9 P 7 9,9 U»4 

32,9 3«. .1 24 9 35. 3 35* '» 

7*5 7-6 3 2 9,% 12*9 

34.4 >5<7 54 5 36 9 3t«6 

7,: 7.1 7 7 5, 6 lt*2 

35, a 37« 1 39,1 30, 4 3e. I 

6 > 6 606 7 , i 8 • 2 9*5 

37c 0 33.5 29-5 39.9 39*6 

6*2 t«. 6 6 7- 6 6*8 

2 6 a 2 39.7 4*. *3 415,3 4l#0 

5 t 7 5,4 4- w 6 9 8.1 

39,2 4 ). 9 42 'I. 42 *5 42a 3 

5.2 5c - 5.-5 6.2 7,5 

4:#2 4.« 9 43 2 43.7 42a5 

4a 8 4,6 4 9 5.6 6 » 6 

4 , #0 4 2o 9 44.2 44*8 44»6 

4o 3 4® 1 4- 3 4,9 5,8 

41, tt *.3,7 45 1 45.6 *5,6 


0.7 0# 7 0.8 1.0 1.2 1.4 1.5 1.7 

34*5 34.1 33,9 33.5 33.0 32.4 31.8 31.2 

0.5 0.5 0.6 0.7 3.8 1*0 1*0 1.1 

34*2 34. C 33.7 33.3 32.8 32.1 31.5 30.9 

0.3 0.3 0.3 0,4 0.4 3.5 0.5 0.6 

34.3 34.0 33.6. 33. 2 32.6 32.3 31.3 3»i*7 

0.0 0.0 0.0 0.0 0.0 0.3 ’3.0 0.0 


. * ce * T * 


3.1 

31.3 

3.4 

30.7 

3.6 

3u.3 

3-. 8 
30c 2 

4.*) 

30.5 

4., 1 

31.4 

2.5 
30. 9 

2.7 

30.4 

2.9 

30.0 

3*1 

30-0 

3.2 

J0.4 

3.) 

31*2 

1.8 

30.6 

2.0 

30. i 

2.2 

29.6 

2. 3 
29,8 

2.4 

30.3 

2,5 

31,2 

1.2 

31.) 

1.3 

29.9 

1 0 5 
29.6 

1*6 
29.. 7 

1.5 

20.2 

1.7 
31 1 1 

0.6 

3C.2 

0.7 

29.7 

0.7 

29.5 

C. 6 
29. 6 

3.8 

30.1 

0.0 

3 l.C 

C.l 

3V..0 

C.O 

29-6 

0*0 

29.) 

0..0 

29.3 

0.3 

29.8 

3. C 
3O.0 

0.32 

3.36 

0.40 

0.44 

0.40 

0.52 


5.6 

32.) 

5,4 

33,7 

5i l 
35,6 

6,7 

37,8 

4® 3 
4V,l 

3. 8 
42,5 

},t> 

44.6 

3 8 
46* C 

4c 3 
46:7 

5o 1 

46. 5 

4*8 

32.4 

4,7 

33,9 

4.5 

35.9 

4. 1 
38.2 

3c 7 
40,6 

3,3 
43c 1 

3.1 
45c 2 

3*2 
46- 7 

3:7 

47.5 

4,4 

47.3 

4» ; 
32.5 

4,7 

34,1 

3c 6 
36,1 

3.5 

3e.5 

3,2 
4*. i 

2*8 

43,6 

2* 6 
45nT 

2 6 
47.4 

3.* 

48,2 

2.7 

49.1 

3*3 

32.6 

3,2 

34,2 

3,1 

36,4 

2.6 

36,8 

:.6 
4iii 4 

2r 3 

44,7 

2. * 
44,2 

? 1 
47.9 

2,4 
49 .£ 

3.0 

48-7 

2c 5 
32.5 

7,4 

34,3 

2.3 

36.5 

2.1 

39,0 

1.9 

41.7 

1.7 

44.3 

1.5 

46,6 

5 

48,4 

1,6 
49 ) 

2,3 

49.2 

1.7 

32.5 

i*5 

34,3 

1.5 
)&, 6 

1.4 

39,2 

1.3 

41.9 

1.1 

44. 6 

..0 

46.9 

TcC 
40 •« 

1.2 
49, 7 

1.6 

49.5 

• c 6 
32.4 

3,8 

34,3 

0.8 

56.6 

0,7 

39,2 

0,6 

42.0 

*' ? 6 

44,6 

jo 5 
47-2 

'-•5 

49,1 

0-.6 

50.0 

0.9 

49.8 

32! 2 

0,0 

34,1 

0.5 

36.5 

f.O 

39.3 

0.0 

42.2 

0.0 

45.0 

0. J 

47.6 

0.0 

46. 6 

0, 0 
5). 4 

0,0 

50,0 

C. 56 

0.60 

0,64 

3.68 

0,72 

0.76 

0.80 

:»84 

?,ce 

0.9} 


EtCH etOCK IS IUDUI 01Si»LiCcHCNT in MICRO niH-'.S . 

2 013*»lACfcM£»4T in HICROINCM .S 

) 
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